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Abstract. We consider a strictly hyperbolic, genuinely nonlinear system of conservation laws 
in one space dimension. A sharp decay estimate is proved for the positive waves in an entropy 
weak solution. The result is stated in terms of a partial ordering among positive measures, using 
symmetric rearrangements and a comparison with a solution of Burgers' equation with impulsive 
sources. 



1 - Introduction 

Consider a strictly hyperbolic system of n conservation laws 



and assume that all characteristic fields are genuinely nonlinear. Call X\(u) < ■ ■ ■ < X n (u) the 
eigenvalues of the Jacobian matrix A(u) = Df(u). We shall use bases of left and right eigenvectors 
li(u), Ti{u) normalized so that 



Given a function u : JR i— > M n with small total variation, following [BC], [B] one can define the 
measures [i % of i-waves in u as follows. Since u 6 BV, its distributional derivative D x u is a Radon 
measure. We define \x l as the measure such that 



ut + f{u) x = 



(1.1) 





if i = j, 
if i^j. 



(1.2) 



[i l = k{u) ■ D x u 



(1.3) 



restricted to the set where u is continuous, while, at each point x where u has a jump, we define 

S({x})=a l , (1.4) 



where Uj is the strength of the z-wave in the solution of the Riemann problem with data u~ = u(x— ), 
u + = u(x+). In accordance with (1.2), if the solution of the Riemann problem contains the 
intermediate states u~ = Uq,U\, . . . ,u n = u + , the strength of the i-wave is defined as 

<Ti = - Xi(ui-i). (1.5) 

Observing that 

a, = k(u + ) ■ (u + - u~) + 0(1) • \u + - u~\ 2 , 
we can find a vector U(x) such that 

\k(x) -k(u(x+)) \ = (D{1) ■ \u(x+) - u(x-)\, (1.6) 

<7j = k(x) ■ (u(x+) - u(x-)). (1.7) 
We can thus define the measure [x 1 equivalently as 

^ = k-D x u, (1.8) 

where h(x) = /^(n(x)) at points where u is continuous, while U(x) is some vector which satisfies 
(1.6)-(1.7) at points of jump. For all x G M there holds 

k(x) - k(u(x))\ =C?(l)-|ti(a:+)-«(x-)|. (1.9) 

We call /U I_ respectively the positive and negative parts of so that 

^ = fi i+ -^-, 1^1 = ^+ + ^-. (1.10) 

It is our purpose to prove a sharp estimate on the decay of the density of the measures fj, l+ . This 
will be achieved by introducing a partial ordering within the family of positive Radon measures. 
In the following, meas^) denotes the Lebesgue measure of a set A. 

Definition 1. Let fi, fj! be two positive Radon measures. We say that fj, ^ // if and only if 

sup < sup /J-'(B) for every s > . (1-H) 

meas(A)<s meas(B)<s 

In some sense, the above relation means that // is more singular than [i. Namely, it has 
a greater total mass, concentrated on regions with higher density. Notice that the usual order 
relation 

H < n' if and only if n(A) < n'(A) for every A C M (1.12) 
is much stronger. Of course \i < \i! implies \i < /i', but the converse does not hold. 
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Following [BC], [B], together with the measures /x* we define the Glimm functionals 



n«)=zyiw, (i-i3) 



Q(u) = ^2(\^\<»\^\){(x,y); x < y) + £ ® H){(x,y) ; x ^ y} . (1.14) 

i<j i 

Let now u = u(t,x) be an entropy weak solution of (1.1). If the total variation of u is small and 
the constant Co is large enough, it is well known that the quantities 

Q(t) = Q(u(t)) , T(t) = V(u(t))+C Q(u(t)) (1.15) 

are non-increasing in time. The decrease in Q controls the amount of interaction, while the decrease 
in T controls both the interaction and the cancellation in the solution. 

An accurate estimate on the measure of positive i-waves in u(t, ■) will be obtained by a 
comparison with a solution of Burgers' equation with source terms. 

Theorem 1. For some constant n and for every small BV solution u = u(t,x) of the system (1.1) 
the following holds. Let w = w(t, x) be the solution of the scalar Cauchy problem with impulsive 
source term 

w t + (w 2 /2) x = -Ksgn(x) ■ j t Q{<t)) , (1.16) 
w(0,x) = sgn(x) • sup ^° ^ . (1.17) 

meas(A)<2\x\ * 

Then, for every t > 0, 

ti + <D x w{t)- (1-18) 



As shown in the next section, the initial data in (1.17) represents the odd rearrangement of 
the function Vi(x) = /i + (] — oo,rc]). The above theorem improves the earlier estimate derived 
in [BC]. For a scalar conservation law with strictly convex flux, a classical decay estimate was 
proved by Oleinik [O]. In the case of genuinely nonlinear systems, results related to the decay of 
nonlinear waves were also obtained in [GL], [LI], [L2], [BG]. An application of the present analysis 
will appear in [BY], where Theorem 1 is used to estimate the rate of convergence of vanishing 
viscosity approximations. 



2 - Lower semicontinuity 

3 



Let ^bea positive Radon measure on M, so that n = D x v is the distributional derivative of 
some bounded, non-decreasing function v : 1R h- ► M. We can decompose 



as the sum of a singular and an absolutely continuous part, w.r.t. Lebesgue measure. The absolutely 
continuous part corresponds to the usual derivative z = v x , which is a non- negative L 1 function 
defined at a.e. point. We shall denote by z the symmetric rearrangement of z, i.e. the unique even 
function such that 

z(x) = z(-x), z{x)>z{x) if 0<x<x', (2.1) 

meas ^{x ; z(x) > c}^ = meas ^{x ; z(x) > c}^ for every c > . (2.2) 

Moreover, we define the odd rearrangement of v as the unique function v such that (fig. 1) 

v(-x) = -v(x) , v(0+) = ^ sing (iR) , (2.3) 

v( x ) = u(0+) + f z(y) dy for x > . (2.4) 
Jo 

By construction, the function v is convex for x < and concave for x > 0. 




The relation between the odd rearrangement v and the partial ordering (1.10) is clarified by 
the following result, which is an easy consequence of the definitions. 

Proposition 1. Let fi = D x v and fj,' = D x v' be positive Radon measures. Call v,v' the odd 
rearrangements ofv,v', respectively. Then fi ^ D x v ■< fi and moreover 

v(x) = sgn(x) • sup , (2.5) 

meas(A)<2\x\ 1 



M ^ A*' ora/y «/ < v'(a;) /or all x > . 



(2.6) 



Two more results will be used in the sequel. By the restriction of a measure fj, to a set J, we 
mean the measure 

0^)04) = M^n J). 

Proposition 2. Let fi,fi' be positive measures. Consider any finite partition M = J\ U • • • U Jat. 
7/i/ie restrictions of /i,// to eac/i set satisfy \i\Jt ■< l^'[Je, then 

Proposition 3. Assume that ji ■< D s w for some nondecreasing odd function w. If — /x|(iR) < e, 

r <s 
/J diD s w + sgn(s) • - 



The next result is concerned with the lower semicontinuity of the partial ordering ^ w.r.t. weak 
convergence of measures. 

Proposition 4. Consider a sequence of measures \i v converging weakly to a measure /j,. Assume 
that the positive parts satisfy //+ ^ Dw v for some odd, nondecreasing functions s i— ► w v (s), concave 
for s > 0. Let w be the odd funtion such that 

w(s) = liminf w u (s) for s > . 

I/-»CXD 

Then the positive part of fi satisfies 

H + <D s w. (2.7) 

Proof. By possibly taking a subsequence, we can assume that w v (s) — > w(s) for all s / 0. 
Moreover, we can assume the weak convergence 

for some positive measures /2~. We thus have 

fj, = fi + — fi~ , n + < £l + , /U~ < /U~ • (2-8) 

By (2.8) it suffices to prove that /i + ^ D s f«, i.e. 

meas (A) < 2s =^> A + (A) < 2iw(s) , (2.9) 
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for every s > and every Borel measurable set A C M. If (2.9) fails, there exists s > and a set 
A such that 

meas(^) = 2s, P + (A) > 2w(s) = 2 lim w v (s). 

v— s-oo 

Since w is continuous for s > 0, we can choose an open set A' ^ A such that, setting s' = 
meas(A')/2, one has 2w(s') < jl + (A). By the weak convergence /i+ — 1 fi + one obtains 

/2 + (A') < Iiminf/i+(A') < 2w(s') < 

reaching a contradiction. Hence (2.9) must hold. □ 

Toward the proof of Theorem 1 we shall need a lower semicontinuity property for wave mea- 
sures, similar to what proved in [BaB]. In the following, Co is the same constant as in (1.15). 

Lemma 1. Consider a sequence of functions u v with uniformly small total variation and call 
the corresponding measures of positive i-waves. Let s h- ► w v (s), v > 1, be a sequence of odd, 
nondecreasing functions, concave for s > 0, such that 



w u + C sgn(s)(Q -Q(u u )) (2.10) 



for some Qq. Assume that u v — ► u and w v — > w in Lj^. Then the measure of positive i-waves in 
u satisfies 



i+ -< D 



w + C oS gn(s)(Q -Q(u)) . (2.11) 



Proof. The main steps follow the proof of Theorem 10.1 in [B]. 

1. By possibly taking a subsequence we can assume that u v (x) — ► u(x) for every x and that the 
measures of total variation converge weakly, say 

= \D x u„\ -^fjfi (2.12) 

for some positive Radon measure /A In this case one has /i" > in the sense of (1.12). 

2. Let any e > be given. Since the total mass of ^ is finite, one can select finitely many points 
Ui, - ■ ■ ,Vn such that 

^({x})<£, for all x£{y u ...,y N }. (2.13) 

We now choose disjoint open intervals Ik = }yk — p, Uk + p[ such that 

^(I k \{y k })<^ k = l,...,N. (2.14) 



Moreover, we choose R > such that 

N 



|J I k C [~R,R], - oo, -i?] U [i?, oo [) < e. (2.15) 

fc=i 

Because of (2.13), we can now choose points po < —R < Pi < ■ ■ ■ < R < p r which are continuity 
points for u and for every u u , such that 

M"({Pfc})=0, u v (p h )^u{p h ) forall h = 0,...,r, (2.16) 
and such that either 

p h -p h -i<—, Ph-i<Vk<Ph, \ph-i,Ph] C 4, (2.17) 
for some fc € {1, . . . , iV}, or else 

lA*l([Ph-i,Pfc]) < A**([Ph-i,Ph]) < £■ (2-18) 
Call Jh = [f^-ijPh]- If (2.18) holds, by weak convergence for some v sufficiently large one has 

\n„\(J h )<£ for all v > vq. (2.19) 

On the other hand, if (2.17) holds, from (2.14) it follows 

H ( J h \ {y k }) < m" (Jh \{y k })<^. (2.20) 

In the remainder of the proof, the main strategy is as follows. 

• On the intervals Jh(k) containing a point yk of large oscillation, we first replace each u v by 
a piecewise constant function u v having a single jump at yk- The relations between the 
corresponding measures \i l v and fil are given by Lemma 10.2 in [B]. Then we take the limit as 
v — > oo. 

• On the remaining intervals Jh with small oscillation, we replace the left eigenvectors li(u u ) by 
a constant vector li(u* h ). Then we use Proposition 4 to estimate the limit as v — > oo. 

3. We first take care of the intervals Jh containing a point yk of large oscillation, so that (2.17) 
holds. For each k = 1, . . . , TV, let h = h(k) G {1, . . . , r} be the index such that yt € Jh = [Ph-i, P/i]- 
For every z/ > 1 consider the function 

u v (x) if x £ U k J h ( k ), 

u v (x) = { u v (ph( k )-\) if x e]Ph( k )-i, Vk[, 
u v {ph) if x € [j/ fc , Ph(fc)]. 



Observe that all functions u, u v are continuous at every point p ,...,p r and have jumps at 
yi,...,yN- Call ftl, i = l,...,n, the corresponding measures, defined as in (1.8) with u re- 
placed by u u . Clearly fll = nl outside the intervals Jh(k) of large oscillation. By Lemma 10.2 at 
p.203 in [B], there holds 



Q(u v ) < Q(u u ) 



V(u v ) + C Q(u v ) < V(u v ) + C ■ Q(u„), 



K + {IR)-ti + (lR) < C [Q(u v )-Q(u v ) 
As a consequence, from (2.10) we deduce 

K + r< D s \t"w v + C sgn(s) (Qo - Q{u v )) 



(2.21) 



where 



T e w(s) 



w(s + e/2) 



if s > 0, 



w(s-e/2) if s<0. 



Indeed, all the mass which in lies on the set 



N 



& = U A(fc) , 

fc=l 



Jfc = [Pfc_l, Pfc] 



is replaced in by point masses at yi, . . . ,Vn- We obtain (2.21) by observing that, by (2.17), 
meas(Q) < e. Moreover, the increase in the total mass is < Co \Q(u v ) — Q(u v )\ . 
Since u v {ph) — »• u{ph) for every h, there holds 

V l {{Vk}) - &{{Vk}) = 0(1) ■ {\u(y k -) - «(p/,(fc)_i)| + \u(y k +) - u(p h(k) )\ 

+ - Uu(Ph{k)-i)\ + \u(p h {k)) ~ Uv(j>h(k))\} (2.22) 



0(1). 



N 



for each k = 1, . . . , N and all v sufficiently large. By construction we also have 



lALK^(fc) \ ivk}) = °> \A{J h {k) \ {y k }) = oil) ■ 



N 



(2.23) 



4. Next, call S = {h; //"(Jj^) < e} the family of intervals where the oscillation of every u v is 
small, so that (2.18) holds. If h G <S, for every x, y € Jh and j/ sufficiently large, one has 

\u v (x) - u u {y)\ < \Hv\{Jh) < e, 

< |/i|(J h )</i lt (Jh)<e. 
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Set u* h = u{ph). By the pointwise convergence u v {ph) — > ^(p/i) and the two above estimates it 
follows 

\u„(x) — u* h \ < e, \u(x) — u* h \ < e, for all xG 4. (2.24) 

5. We now introduce the measures such that 

restricted to each interval Jh, h <E S where the oscillation is small, while 

fC = K 

on each interval Jh = Jh(k) where the oscillation is large. Observe that the restriction of jX\, to 
Jh(k) consists of a single mass at the point y^. Namely, fil({yk}) is precisely the size of the i-th 
wave in the solution of the Riemann problem with data u~ = u l/ (ph(k)-i), u + = u v (ph(k))- 
We define w v as the non-decreasing odd function such that 

w v (s) = sup ; s>0 . (2.25) 

meas(A)<2s ^ 

By possibly taking a further subsequence we can assume the convergence 

Q{u v )^>Q, /}*->•/}*, w u (s) -> w(s) . 

Using (2.16), we can apply Proposition 4 on each interval Jh and obtain 

fi i+ <D s w. (2.26) 

6. Observe that, by (2.24) and (2.19), 

| At - m1|(<4) = 0(1) -e^(Jh) heS, (2.27) 
From (2.21) and the definition of w v at (2.25) it thus follows 

w v (s) < T s w v {s) + C [Qo - Q{u v )] + 0(1) ■ e s>0. (2.28) 
Letting i/-»oowe obtain 

w(s) <T £ w(s) + Co[Qo-Q]+0(l)-£ s>0, (2.29) 

Q = lim Q(u„) > lim Q(u u ) - 0(1) • e > Q{u) - 0(1) • e , (2.30) 

V— >CO I/-KX) 

because of the lower semicontinuity of the functional u i— > Q(-u). From (2.26), (2.29) and (2.30) we 
deduce 



li 



i+ -< D 



T £ w + sgn(s) (C [Q - Q{u)} + 0(1) • e 
By (2.22)-(2.24), our construction of the measure fi l achieves the property 

-fi i+ \(M) = 0(1) -e. 



Hence, by Proposition 3, 



T £ ™ + sgn(s) (C [Q - Q(u)] + 0(1) • e) 



Since e > was arbitrary, this proves (2.11). □ 
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3 - A decay estimate 



The second basic ingredient in the proof is the following lemma, which refines the estimate in 
[BC]. 

Lemma 2. For some constant k > the following holds. Let u = u(t,x) be any entropy weak 
solution of (1.1), with initial data u(0,x) = u(x) having small total variation. Then the measure 
/4 + of positive i-waves in u(t,-) can be estimated as follows. 
Let w : [0, r[ xlR h- > M be the solution of Burgers' equation 

w t + (w 2 /2) x = (3.1) 

with initial data 

w(0,x) = sgn(x) • sup ^ (A) . (3.2) 

meas(A)<2\x\ 2 

Set 

w(t, x) = w(t- , x) + /«sgn(x) • [Q(u) - Q(u(t)) . (3.3) 

Then 

& < D x w{t) . (3.4) 



Proof. The main steps follow the proof of Theorem 10.3 in [B]. We first prove the estimate (3.3) 
under the additional hypothesis: 

(H) There exist points yi < ■ ■ ■ < y m such that the initial data u is smooth outside such points, 
constant for x < y\ and x > y m , and the derivative component li(u)u x is constant on each 
interval ]yg, Moreover, the Glimm functional 1 1— ► Q(u(t)) is continuous at t = r. 

1. The solution u = u(t, x) can be obtained as limit of front tracking approximations. In particular, 
we can consider a particular converging sequence (u v ) v >i of e„-approximate solutions with the 
following additional properties: 

(i) Each i-rarefaction front x a travels with the characteristic speed of the state on the right: 

x a = Xi(u(x a +)). 

(ii) Each i-shock front x a travels with a speed strictly contained between the right and the left 
characteristic speeds: 

Xi(u(x a +)) < x a < Xi(u(x a -)). (3.5) 
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(iii) As v — > oo, the interaction potentials satisfy 

Q(u„(0,-)) -<3(«). (3.6) 

2. Let -u^ be an approximate solution constructed by the front tracking algorithm. By a (general- 
ized) i- characteristic we mean an absolutely continuous curve x = x(t) such that 

x(t) G [Xi(u v (t,x—)), \i(u„(t,x+))] 

for a.e. t. If u v satisfies the above properties (i)-(h), then the i-characteristics are precisely the 
polygonal lines x : [0, r] i— > M for which the following holds. For a suitable partition = to < h < 
■ ■ ■ < t m = t, on each subinterval [tj-i, tj] either x(t) = Xi(u v (t,x)), or else x coincides with a 
wave-front of the i-th family. For a given terminal point x we shall consider the minimal backward 
i- characteristic through x, defined as 

y(t) = min {x(t) ; x is an z-characteristic, x(t)=x}. 

Observe that y(-) is itself an ^-characteristic. By (3.5), it cannot coincide with an i-shock front of 
u on any nontrivial time interval. 

In connection with the exact solution u, we define an ^-characteristic as a curve 

t h- ► x(t) = lim x v (t) 

which is the limit of i-characteristics in a sequence of front tracking solutions u v — > u. 

3. Let e > be given. If the assumption (H) holds, the measure ^ of i-waves in u[r) is supported 
on a bounded interval and is absolutely continuous w.r.t. Lebesgue measure. We can thus find a 
piecewise constant function -0 r with jumps at points x\{f) <x^{t) < ... < x~n{t) such that 

da; L . c-i V 



dx < e , / - V T ) dx = j = 1, . . . , N - 1 . (3.7) 



X 














r 


2 /\ V*) / \ 






1 X 1 


1 / Vi Ai \ 







figure 2 
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To prove the lemma in this special case, relying on Proposition 2, it thus suffices to find 
i-characteristics t Xj (t) such that the following holds (fig. 2) 

(i) For each j = 1, . . . , N, the function i[) T is constant on the interval ]xj(r), Xj+i(r) [ and (3.7) 
holds. Moreover, either Xj(0) = Xj+i(0), or else the derivative component if) = l % (u)u x (0, ■) 
is constant on the interval ]xj(0), x J+ i(0)[. 

(ii) An estimate corresponding to (3.3)-(3.4) holds restricted to each subinterval [xj(t), xj+i(t) [. 
We need to explain in more detail this last statement. Define 

Ij(t)= [ Xj (t), x j+1 (t)[, A j = {(t,x); t€[0,r], x€lj(t)}. 

For each j, we denote by Tj the total amount of wave interaction within the domain Aj. This is 
defined as in [B], first for a sequence of front tracking approximations u u , then taking a limit as 
v — ► oo. Furthermore, we define the constant values 

Vj=V{x) x€Ij(t), 

^ = ^°(x) x€lj(0), 

Call 

the initial amount of positive i- waves inside the interval Ij . 

For each interval Ij, we consider on one hand the function Wj corresponding to (3.2)-(3.3), 
namely 

wj(a) = min |a° , r + ( ^ 0) _ 1 | + nTj ■ sgn(s) . 

Here (^°) _1 = in the case where Xj(0) = x J+ i(0). This may happen when the initial data has 
a jump at Xj(0), and the corresponding measure has a Dirac mass (with infinite density) at 
that point. 

On the other hand, we look at the nondecr easing, odd function rjj such that 
Tjj (s) = min | s, [xj+i (r) - Xj (r)] } s > . 

Our basic goal is to prove that (fig. 3) 

Vj(s) < for all s > . (3.8) 

Indeed, by (3.7), for s > one has 

sup < r)j{s) + Sj 

meas(A)<2s z 
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with 



Proving (3.8) for each j will thus imply 

±w(t,x) =w{t-, x) + Ksgn(x)- [Q(u) - Q(u(t)) +0(1) -e . 
Since e > was arbitrary, this establishes the lemma under the additional assumptions (H). 



a J +Cr j 



o 



Wj(8) 



Tlj(8) 



figure 3 



4. We now work toward a proof of (3.8), in three cases. 
Case 1: = 0. 

Case 2: Xj(0) = x j+1 (0) and o-° > 0. 

Case 3: Xj(0) < x j+1 (0) and a<] = (x j+1 (0) - Xj(0)) > 0. 

In Case 1 the proof is easy. Indeed, the total amount of positive i-waves in Ij(r) is here 
bounded by a constant times the total amount of interaction taking place inside the domain Aj, 
i.e. 

^ T + {l j {r))<C -Y j 
for some constant Cq. On the other hand 

Wj(s) = kTj ■ sgn(s) . 

Choosing k > Co we achieve (3.8). 
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5. Since Case 2 can be obtained from Case 3 in the limit as Xj+i — Xj — > 0, we shall only give a 
proof for Case 3. 

We can again distinguish two cases. If the amount of interaction Tj is large compared with 
the initial amount of i- waves, say 

then the bound (3.8) is readily achieved choosing n > 8Cq- Indeed, for s > we have 

Vj (s) < \ ^ {Ij (r)) < Co^ + a Q 5 < 7C Tj . 
The more difficult case to analyse is when Tj is small, say 

r 3 -<o-5/6C . (3.9) 

Looking at figure 3, it clearly suffices to prove (3.8) for the single value 

Sj + i(r) -gj(r) 
s ~ s i ~ 2 

Equivalently, calling 

Zj{t) = x j+1 (t) - Xj(t) 

the length of the interval Ij(t) and 

^=^ + {lj(r))=Zj(r)^ 
the total amount of positive i-waves inside Ij(r), we need to show that 

< Tj<2/,r j +min|a J , T+ ^o } -i } • ( 3 - 10 ) 
By the approximate conservation of i- waves over the region Aj , we can write 

<Jj < o-J + C Tj . (3.11) 
Using (3.11) in (3.10), our task is reduced to showing that 

°i ^ 2 ^ + r + ^0)-! ( 3 - 12 ) 
for a suitably large constant k. Because of (3.11), it suffices to show that 

Zj(T)>(0°-C'Tj)(T+(^)-l) 

= [zj(0) + ra° j ]-C'(r + (^)- i )Tj 
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figure 4 



for a suitable constant C. 

6. We now prove (3.13). Notice that, by genuine nonlinearity and the normalization (1.2), if no 
other waves were present in the region Aj we would have Tj = and 

In this case, the equality would hold in (3.13). 

To handle the general case, we represent the solution tiasa limit of front tracking approxi- 
mations Uj,, where for each v > 1 the function u„(0, •) contains exactly v rarefaction fronts equally 
spaced along the interval lj(0). Each of these fronts has initial strength <r Q (0) = o^/v. For 
a = l,...,v, let y a (t) G Ij(t) be the location of one of these fronts at time t G [0, r], and let 
a a (t) > be its strength. Moreover, call 

J a (t) = [y a (t), y a +i(t)[, A a = {(t,x); t€[0,r], xeJ a (t)}, 

and let T a be the total amount of interaction in u v taking place inside the domain A a . 
We define a subset of indices X C {1, . . . , u} by setting a G X if 

5C r a ><7 a (0) = o-J/"- (3.14) 

Observe that, if a £ T, then 

cr a (t) 



<r Q (0) 



< - for all t G [0, r] 
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In particular, if a, a + 1 ^ 1, then the interval J a {t) is well defined for all t € [0, r]. Its length 



satisfies the differential inequality 

d_ 

dt 



z a {t) = y a +i(t) - y a (t) 



^a(t)>W a (t)-Ci- Yl M ( 3 - 15 ) 



/3ec«(t) 

for some constant Ci. Here 

W a (t) = [amount of i- waves inside the interval J a (t)] 
>a a (0)~C T a , 



(3.16) 



while C a (t) refers to the set of all wave fronts of different families which are crossing the interval 
J a at time t. Calling W' a the total amount of waves of families ^ i which lie inside J a (0), we can 
now write 

f ( E k/9l]*<(!j«^ (3.i7) 

Indeed, by strict hyperbolicity, every front ap of a different family can spend at most a time 
0(1) • -Zq inside J a . Either it is located inside J a already at time t = 0, or else, when it enters, 
it crosses y a or y a +\. In this case, since a, a + 1 ^ X, by (3.14) it will produce an interaction of 
magnitude Ict^UqI > \ap ■ a®\/2v. The second term on the right hand side of (3.17) takes care of 
the new wave fronts which are generated through interactions inside J a . The last term takes into 
account wave front of different families that initially lie already inside J a at time t = 0. Integrating 
(3.15) over the time interval [0,r] and using (3.16)-(3.17) one obtains 

z a (r) > z a (0)+T^-O(l)-TF a -O(l)-(m^zUt))~-r a -O(l)-( Zji0) + 1 ) W' a . (3.18) 
v \te[o,r] J a- \ v J 



7. To proceed in our analysis, we now show that 

max z a (t) < 2z a (r) . (3.19) 
te[o,r] 

Indeed, let r' G [0, r] be the time where the maximum is attained. If our claim (3.19) does not 
hold, there would exists a first time t" E [t',t] such that z a {r") = z a {r')/2. ^From (3.15) and 
the assumption W a (t) > it follows 

f & E W?\dt> Z -^. (3.20) 

pec a (t) 
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Using the smallness of the total variation, a contradiction is now obtained as follows. Call 

$(t) = C Q(t) + £ K{t,x p (t)) \a p \ , 

where the sum ranges over all fronts of strength ap located at x@, of a family kp ^ i. The weight 
functions <j>j are defined as 



<f)j(t,x) = < 





y a +i{t) - x 

y a +i(t) - y a (t) 



if x>y a+1 (t), 

if x e [y a (t), y a +i(t)] , 

if x<y a (t), 



in the case j > i, while 

4>j(t,x) 



x - y a (t) 



if x>y a+ i(t), 
if x € [y a (t), y a +i(t)] , 
if x<y a (t), 



y a +i(t) - y a (t) 


in the case j < i. Because of the term CoQ(t), the functional $ is non-increasing at times of 
interactions. Moreover, outside interaction times a computation entirely similar to the one at 
p. 213 of [B] now yields 



dt 



z(t) 



(3.21) 

jii i 

PeC a (t) 

for some small constant c > related to the gap between different characteristic speeds. From 
(3.20) and (3.21) respectively we now deduce 



/3ec a (t) 



// / 



/3ec„(t) 



d$(t) 



c 



CO 



Since <&(i) = 0(1) ■ Tot.Var.{u(i)}, by the smallness of the total variation we can assume 3>(t') < 
2C 1 /c . In this case, the two above inequalities yield a contradiction. 

8. Using (3.19), from (3.18) we obtain 
z i( r ) = E Zq ( t ) - E z «( r ) 



Ka<u 



a<£2 



- E ] i + c 2 (^°)r Q 



0(l)-rr,. -0(1) 



(3.22) 



CTn ^ -i). Trj - (i).fi!fi±l 
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By (3.14) the cardinality of the set X satisfies 



hence 



In turn, this implies 



— n i 

°3 



£ U(0) + r^) > (^(0)+ra°) (l - ^ > (z^+ra*) -SCoY^Yj-SCotYj . (3.23) 



Using (3.23) in (3.22), observing that 



Zj (0) _ x j+1 (0) - Xj (0) 



and letting i/^oowe conclude 



Zj(r) > { Zj (0) + r<x°) - 0(1) • ft®" 1 ^- - 0(1) ■ rY . 
This establishes (3.13),for a suitable constant C . 

9. In the general case, without the assumptions (H), the lemma is proved by an approximation 
argument. We construct a convergent sequence of initial data u v — > u which satisfy (H) and such 
that 

u u ^u, Q(u u ) -» Q(u) , |f**+ - ^ + | -» 0. 
Calling the solution of (3.1) with initial data 

1^(0,2) = sgn(x) • sup , 

mea.s(A)<2\x\ * 

by the previous analysis we have 



H% v ^ D x w v {t v -) + sgn(x) • [Q{u v ) - Q(u v {t v ))] 

Observe that w v (r v —) — > w(t—) in Lj oc . Choosing /« > Co, by the lower semicontinuity result 
stated in Lemma 1 we now conclude 



/4 + < D a 



w(t-) + Ksgn(x) • [Q(u) - Q(u(t))] 



□ 
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4 - Proof of the main theorem 



Using the previous lemmas, we now give a proof of Theorem 1. For a given interval [0,r], 
the solution of the impulsive Cauchy problem (1.17)-(1.18) can be obtaines as follows. Consider 
a partition = t < ti < ■ ■ ■ < = t. Construct an approximate solution by requiring that 

w(0,x) = Vi(x), 

w t + (w 2 /2) x = (4.1) 

on each subinterval [t k -i,t k [, while 

w(t k ,x) = w(t k - , x) + Ksgn(x) • [Q(ifc_i) - Q(t k )] • (4.2) 

We then consider a sequence of partitions = t% < t\ < ■ ■ ■ < t v Nv = r, and the corresponding 
solutions w v . If the mesh of the partitions approaches zero, i.e. 

hm sup | = 0, 

i/— oo k 

then the approximate solutions w v converge to a unique limit, which yields the solution of (1.17)- 
(1.18). 

Call T the set of nondecreasing odd functions, concave for x > 0. This set is positively 
invariant for the flow of Burgers' equation (4.1). Moreover, this flow is order preserving. Namely, 
if w,w' G T are solutions of (4.1) with initial data such that w(0,x) < w'(0,x) for all x > 0, then 
also 

w(t,x) < w'(t,x) for all t, x>0. 

Equivalently, 

D x w(0) r< D x w'(0) D x w{t) ± D x w'{t) 

for every t > 0. For each fixed v, we can apply Lemma 2 on each subinterval [t^_ 1 ,t k r \ and obtain 

By induction on fe, this yields 

/4 + :< AcMr), ( 4 - 3 ) 

where is the approximate solution constructed according to (4.1)-(4.2). Letting v — > oo and 
using Lemma 1 , we achieve a proof of Theorem 1 . □ 
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